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We put forward a reconstruction scheme for dark energy which employs a number of observation- 
ally motivated geometrical constraints expected to be satisfied by the luminosity-distance relation. 
This generalises a scheme recently put forward, with the important additional feature that it allows 
a wider range of possible equations of state to be explored. 

Applying this scheme to the recent SNLS supernovae data we find a number of interesting results: 
(i) The ACDM model is compatible with the data within the 2a significance level, but not at 
the la level, (ii) Allowing the dark energy model to have an equation of state which is free to 
become phantom in the redshift interval [0, 1], we find a sharp increase in the degeneracy in the 
reconstructed equation of state for redshifts greater than ~ 0.6, brought about by a sudden widening 
of the confidence contours. We trace the origin of this enhancement in degeneracy to the fact that 
in such models the equation of state can diverge. This is important in order to devise strategies 
to reduce this degeneracy, specially in the light of near future observations. There is, however, no 
such large degeneracies for dark energy models whose equations of state remain larger than —1 in 
the redshift interval [0, 1]. (iii) There exist reconstructed solutions which fit the data well and show 
transient acceleration, with a very recent phase of deceleration. 

PACS numbers: 04.50.-|-h, 98.80.Es 



I. INTRODUCTION 

Recent observations have provided strong evidence to 
suggest that the universe has recently undergone a phase 
of accelerated expansion. These include evidence from 
supernovae data [l|, 0] as well as CMB measurements 
d, 01 and observations of large scale structure [1, [1] . 

The simplest way to explain this late acceleration is 
by invoking the presence of a positive cosmological con- 
stant. Although this is compatible with all observational 
constraints to date, it is not clear how to reconcile the 
presence of such a small parameter with fundamental the- 
ories. As a result a range of other scenarios have been 
considered to account for this acceleration. An important 
class of these scenarios - the so called 'dark energy' sce- 
narios - models the underlying accelerating mechanism in 
terms of an effective perfect fluid with a negative equa- 
tion of state (EOS). The attractive feature of these mod- 
els is that they can, in addition to cosmological constant, 
mimic a wide range of scenarios with a variable EOS, 
including quintessence models, braneworld models as 
well as models based on loop quantum cosmology . 

An urgent question at present is how to distinguish be- 
tween these possibilities. Two approaches are possible in 
principle. One may choose a particular model and com- 
pare its predictions with observations. Alternatively, one 
can try to reconstruct the properties of the 'dark energy'. 
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including its equation of state, directly from observations. 

A number of attempts have recently been made at such 
reconstructions, employing different data sets, including 
supernovae data as well as CMB and baryon acoustic 
oscillation (BAO) data (see [l^, [Hi for a recent review 
and references therein). Among these schemes are the 
so called model-independent ones which attempt to re- 
construct the cosmological parameters directly from the 
data without makin g p articular assumptions regarding 
their forms [13, [3, Il5| . It is important to note that 
the interpretation of observational data always requires 
a theoretical framework. In that sense all reconstructions 
are framework dependent. The question is how wide a 
set of models such a framework allows. 

The aim of this paper is to generalise the reconstruc- 
tion scheme recently put forward in [^ by allowing a 
wider set of equations of state to be explored. 

For our reconstructions we shall employ supernova 
data recently released by the Supernova Legacy Survey 
(SNLS) [2,' as well as the baryon acoustic oscillation 
data (BAO) [l^. Using this scheme we shall reconstruct 
the luminosity-distance curves best fitting the data and 
hence the corresponding cosmological parameters includ- 
ing the Hubble parameter and the EOS. We find that the 
generalised scheme leads to a number of interesting re- 
sults (see below). 

The plan of the paper is as follows. In Section [ll] we 
give a brief description of the reconstruction scheme and 
the data sets used in the reconstructions. In Section Hill 
we summarise the results of our reconstructions. This 
includes a brief summary of our reconstruction of the 
mock data in subsection IIII Al a summary of our re- 
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construction of the luminosity-distance relation and its 
derivatives as well as the Hubble parameter in subsec- 
tion [IIIbI an estimate for the cold dark matter density 
parameter flmg using the BAO data in subsection IIII B 2l 
and finally a summary of our reconstructions of the EOS 
and the deceleration parameter in subsection IIII B31 In 
Section IIVI we explain the origin of the degeneracy (i.e. 
sudden widening of the confidence contours) appearing in 
the reconstruction of the EOS and discuss how it can be 
reduced. In Section|V|we consider two important subsets 
of reconstructed curves representing respectively models 
with quintessence dark energy and transient acceleration. 
Finally, Section lVTl contains a summary of our results and 
conclusions. 



II. THE RECONSTRUCTION SCHEME 

In our reconstruction scheme the luminosity-distance 
di{z), expressed in terms of the Hubble parameter H and 
redshift z as 

di{z) = c{l + z) -, (1) 

is reconstructed subject to two geometrical constraints 
expected to be satisfied by the luminosity-distance rela- 
tion. These are: 

[A ] rf; > 

[B ] d'l > 0. 

where a prime denotes differentiation with respect to the 
redshift z. 

The assumption A is simply a consequence of the ob- 
served expansion of the universe. Assumption B can be 
shown to imply that the deceleration parameter q < I. 
This would be satisfied by models whose current equa- 
tion of state satisfy wq < 0.46 (see below). These include 
non-accelerating models as well as the important class of 
models that are currently accelerating (with q < 0) and 
which in the past tend to an Einstein de-Sittter model 
(with q — 1/2). We shall show that the assumption q < I 
is preferred by the data at Scr for z < 0.5, whereas at 
larger redshift q could be greater than 1 due to large 
degeneracy if it was not constrained by condition [B]. 

Our reconstruction scheme generalises a scheme re- 
cently considered by the authors [9j , by dropping a third 
geometrical constraint, namely d"' < 0. The motivation 
for this constraint was that it is satisfied by the Einstein- 
de Sitter and ACDM solutions, which constitute impor- 
tant past and future phases of the 'standard' cosmolog- 
ical model. Despite these nice asymptotic features, this 
assumption, though compatible with observations, is not 
directly imposed by them. One can in fact consider evo- 
lution histories for the universe which are compatible 
with current observations and which do not satisfy this 
assumption for all times. For example, this constraint 
is not respected by some equations of state commonly 



considered in the literature, such as w = wq + wiz or 
w — Wq + Wizjil + z), where Wq and wi are real con- 
stants. It is thus interesting to study the consequences 
of dropping this assumption on the results of the recon- 
structions. 

To proceed we shall employ as our theoretical frame- 
work General Relativity with a dark energy component 
treated as a perfect fluid with a density p and an EOS 
w{z). The Hubble function H is in this case given by 

= Hi [a„o(l + zf + nop/po] , 

where Q,m and f2 are the dark matter and dark energy 
density parameters respectively and the subscript " 0" in- 
dicates the values quantities take at the present time. 
We shall employ the generalised scheme discussed above, 
based on constraints [A] and [B] , in order to reconstruct 
the cosmological parameters. The data used in this study 
are the 115 supernovae with redshifts z < 1.01 from 
the Supernova Legacy Survey (SNLS) To constrain 
we also employ the baryon acoustic oscillation peak 
(BAO) recently detected in the correlation function of 
luminous red galaxies (LRG) in the Sloan Digital Sky 
Survey M- 

To test the goodness of fit of our reconstructions 
we shall, in line with many other works in this area 
[Tsl . [TtI . [20| , use the standard minimization. For the 
supernovae data, is defined by 

2 _ A (m°^^- -m*''0^ 



where n is the number of data points, m°'"* and to*'' are 
respectively the observed and the theoretical magnitudes 
(obtained from the reconstructed luminosity-distance re- 
lations) and (7i is the uncertainty in the individual to°^*' . 

III. RECONSTRUCTIONS 

Our reconstruction procedure begins by generating a 
large set of luminosity-distance curves which satisfy the 
constraints [A] and [B] while fitting the SNLS supernovae 
data (see next Section for details). Each luminosity- 
distance curve di is constructed using a discrete set of 
equally spaced points, and satisfying the constraints [A] 
and [B], which when discretised respectively take the 
forms 

rfii+i > di, 
d'ii+2 ~ dii^i > rfii+i — di- 

In practice we subdivided the interval z e [0, 1] into eight 
sub- intervals of different length allowing us to recover the 
ACDM model at la level with the mock data (see section 
IIII A[) Taking a larger (smaller) number of sub-intervals 
leads to over-fitting (under-fitting) the mock data, re- 
sulting in the ACDM model not to be recovered at the 
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1(7 level. We have checked that using different choices of 
sub-intervals (with the same number of sub-intervals) 
do not qualitatively change our results. In the figures 
presented here, second order interpolation was used to 
obtain the reconstructed values between the steps. We 
have checked that using third order interpolation does 
not qualitatively change our results. 

Using the set of luminosity-distance curves thus deter- 
mined, we calculate for each curve 



H 



d',- 



di 



1 



1 -1 



1 + z J 1 + z 



and 



H' = -c{l + z) 



2di + {l + z) [-2d[ + (1 + 



0<] 



{dl - (1 + z)d[) 



'/^2 



(2) 



(3) 



which can then be used to reconstruct the cosmological 
parameters, inclduing the equation of state, without as- 
suming specific ansatzes (see Q for further details.) 



A. Reconstruction using mock data 

To test our generalised reconstruction scheme we first 
employed mock data, consisting of supernovae data as in 
(see also [H]). The data was generated by assum- 
ing the supernovae have the same distribution in red- 
shift and error bars as the real (SNLS) data, but with 
the luminosity-distances replaced by those correspond- 
ing to the ACDM model plus a gaussian noise. The 
corresponding di curves were then obtained through in- 
terpolation. The reconstructed curves were chosen such 
that their per degree of freedom (xdof) ^^^^ than 
1.20. The best fit was found to be 113.15 which is to 
be compared with — 114 for the ACDM model. Fig- 
ures[T]show the reconstructions of dj and d[' together with 
their la, 2a and 3(t confidence contours respectively. As 
can be seen, the ACDM model lies within the \a confi- 
dence contours, thus demonstrating that the reconstruc- 
tion scheme is reliable. 
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Figure 1: Reconstruction of and d'^ using the mock data. The Icr, 
la and 3cr confidence contours are shown as dashed curves. The KCDM 
model (solid curve) with \^ — 114 lies within the Icr confidence contour. 
The d[ and d[' curves were reconstructed using an interpolation with a 
regular redshift interval of 0.07. Similar interpolation will be used in 
subsequent figures. 



2. Using BAO data to constrain f2mO 

The supernovae data - and hence the reconstructed 
luminosity-distance curves obtained here - do not con- 
strain tightly the value of rJmO (see Q for an explana- 
tion about degeneracy on this parameter). To further 
constrain fimo we require additional input from obser- 
vations. Here we use supernovae together with baryon 
oscillation data [1^1 to recalculate the values corre- 
sponding to each luminosity-distance curve. The results 
are depicted in Fig. [4l As can be seen, the best fit value 
for r^mo is 0.28, with the la confidence interval given by 
0.26 < fimo < 0.3. We shall use this best fit value in the 
rest of the paper. 



B. Reconstruction using real data 

1. Reconstruction ofdi, d'l, d" and the Hubble function 

In this subsection we briefly summarise the results of 
our reconstructions, using the real data, for the quanti- 
ties dl, d[, d'l' (depicted in Fig. [2]) as well as the Hub- 
ble function (shown in Fig. [3]). As can be seen, the 
ACDM model is compatible with the data at 2a sig- 
nificance level, but not at la level since the best fit 
X^ = 112.11. Although this is not significant enough 
to rule out the ACDM model, it does highlight a con- 
sequence of dropping the constraint on the third deriva- 
tive of the luminosity-distance (i.e. d'/' < 0) used in 0. 
Keeping this constraint has the consequence that ACDM 
model becomes compatible with the data at Icr level. 



3. Reconstruction of the EOS and the deceleration 
parameter 

Our reconstructions of the EOS are depicted in Figs. 
[51 A number of comments are in order. 

The reconstructed EOS seems to have a maximum at 
z = 0. This is indeed expected and is a consequence of 
the fact that we have assumed c?" > 0. This maximum 
value can be readily calculated in the following way. Re- 
call that at present (z = 0) the equation of state is given 
by (see e.g. [l9]): 

(1 — f^mo) 
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Figure 2: Shown arc the reconstructions of di (with SNLS data), d[ 
and d'l' . The ACDAl model, with — 114, is compatible with the 
data at 2(T confidence level, but not at Icr level. 



We also have 



d'l'iz = 0) 



Ho 



Ho 



Now since di{z = 0) > 0, this implies that Hq < 2Ho 
and thus wq < ^^Jq — 0.46. 

The reconstructions of EOS also show a sudden large 
degeneracy, brought about by a sudden widening of the 
confidence contours, at around a redshift of z ~ 0.6. This 
makes it extremely difficult to deduce the nature of EOS 
beyond this redshift. We shall discuss this issue in the 
next subsection. 

An alternative way of testing the closeness of EOS to 
— 1 is to use the weighted avera ge of the equation of 
state, the so called w-probe (iD) The advantages 



Figure 3: Figure showing the reconstruction of the Hubble function 
from real data. The dashed lines indicate the Icr, 2(T and 3cr confi- 
dence contours. The upper black line shows the Hubble function for 
the AC DAI model. The lower black line indicates the values of H/Hq 
for which the EOS diverges (see section IHI B 3t . 




Figure 4: Figure showing the distribution of a-s a function of QmQ 
obtained by using the SNLS and BAO data. 



of this probe are that it is independent of the value of 
fimo and does not depend upon the second derivative of 
the luminosity-distance di. The calculated values of w 
are depicted in Fig. [6l As can be seen the constraints 
on the w are better than those on w, even for larger val- 
ues of z for which H) remains finite while w diverges. The 
w-probe provides an alternative confirmation of our re- 
sults concerning the reconstruction of the EOS, namely 
that the ACDM model is in agreement with the data at 
2(7 confidence level. 

We also reconstructed the deceleration parameter, 
which is depicted in Fig. [T] An interesting feature of 
these reconstructions is that they show some solutions 
can possess a transient acceleration. For example at la 
significance level, decelerated expansion can occur for 
0< z < 0.13 and for z > 0.38, while for 0.13 < z < 0.38 
the results show an accelerating expansion. We shall re- 
turn to this issue in section fV Bl below. 



IV. NATURE OF EOS DEGENERACY AND 
STRATEGIES TO REDUCE IT 

Given the important consequences of such sudden 
widenining of the confidence contours for the determi- 
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Figure 7: Reconstruction of the deceleration parameter with its 1, 2 
and 3(T confidence contour 
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Figure 5: Figures showing the reconstruction of the EOS. The upper 
panel shows the sudden widening of the confidence contours around 
z ~ 0.6. The lower panel shows the blow up of this curve. 
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(see the top panel of Fig. [5|). Thus as the reconstructed 
Hubble function approaches the curve C, the confidence 
contours widen and degeneracies increase. In particu- 
lar, the EOS blows up when the reconstructed Hubble 
function intersects the curve C. Consequently some very 
large values of the EOS can fit the data beyond the red- 
shift z > 0.6, leading to the observed degeneracy. In fact 
as can be seen from the plot of curve C, some recon- 
structed solutions within the la confidence contour lines 
correspond to diverging EOS. It is worth noting here that 
a diverging effective EOS cannot necessarily be ruled out 
as unphysical, since it could be an effective mathemati- 
cal representation of some physical models such as gener- 
alised gravity models as well as braneworld models. 

Having seen the cause of the sudden widening of the 
confidence contours and the resulting large degeneracies, 
it is important to determine how they could be reduced. 
Clearly to decrease the observed degeneracy in the 
EOS, one needs to prevent C from intersecting the 
reconstructed H/Hq curves, which implies that the data 
will have to be made more precise for redshifts beyond 
z = 0.6. 



Figure 6: The reconstruction of the w-probe w. The dotted curves 
show the confidence contours at 1, 2, and 3(7 levels respectively. 



nation of EOS, it is important to understand their origin 
in order to find strategies to reduce them. 

To do this, recall that the dark energy EOS can be 
written as |14| 



2(l+z) H' 
3 H 



1 



This shows that the EOS diverges when H/Hq — 
•\/17mo(l + z)^- To see the role of this divergence for the 
above degeneracies we have plotted this curve (referred 
to as C) in Fig. [3] As can be seen the curve C inter- 
sects the reconstructed curves H/Hq, fitting the data, 
at around z > 0.6, which leads to the sudden widening 
of the confidence contours and the resulting degeneracies 
observed in the reconstruction of the EOS for z > 0.6 



To show this more clearly we proceeded as follows. We 
removed all reconstructed H curves which intersected the 
C curve, as can be seen in Fig. [8l The corresponding di 
curves are shown in Fig. [51 This shows the extent to 
which the accuracy of observations needs to be enhanced 
in order to prevent the intersections of reconstructed di 
curves with the curve C, up to redshifts z = I. The corre- 
sponding reconstructed EOS is shown in Fig. [TOl As can 
be seen the removal of H curves intersecting the C curve 
has reduced the degeneracy for larger redshifts. The re- 
sult is rather modest, however, and a further reduction of 
degeneracies, specially at higher redshifts would require 
more severe narrowing of the reconstructed di cone which 
in turn implies more accurate observations (smaller error 
bars) as well as less dispersion. 

Another important point to note in this regard is the 
relation between the C curve and the ACDM curve (see 
Fig. [5]). It is easy to see that these curves intersect 
at z = oo. More importantly the distance between the 
two curves decreases as z increases. This implies that to 
avoid the intersections (and the resulting degeneracies) 
the data has to become more and more accurate as z 
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increases, which is of course a very difficult constraint to 
meet observationally. 
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Figure 8: Reconstruction of the Hubble parameter H with tlic re- 
constructed curves intersecting the C curve removed. The upper solid 
curve represents the ACDM curve and the lower one the C curve. 
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Figure 9: Reconstruction of the di curves corresponding to the H 
curves depicted in Fig. [s] 
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Figure 10: Reconstruction of the EOS corresponding to the recon- 
structed curves in Figs. [8]and|9] 



for details). The reconstructed kinetic and potential en- 
ergies of the scalar field, with their corresponding con- 
fidence contours, are depicted in Figs. [Tl] and [121 As 
can be seen from these reconstructions, the kinetic and 
potential energies are well constrained over the interme- 
diate times in the range z G [0, 1] while they are poorly 
constrained at present (z ^ 0) and at early times (z ~ 1). 
The former is due to the fact that we have no data in fu- 
ture and the latter is the result of the fact that we have 
not used any data beyond the redshift of z > 1 and that 
the error bars are larger for earlier data. 

The best value found is 114, i.e. very close to that 
for the ACDM model, indicating that we do not find 
models that fit the data qualitatively better than the 
ACDM model. Thus ACDM model belongs to the la 
confidence contour. 

We note that the quintessence dark energy does not 
suffer from the degeneracy discussed in connection with 
the dark energy EOS in section UlI B 31 To see why, con- 
sider two functions of z, /i and /2, say. If /( > /2 for 
all z and initially /i(0) > /2(0), then /i > /2. Now, 
the dark energy density p can in general be written as 
p/po = J3{w{z)+i)/{i+z)dz^ ^Yiih for the ACDM model 

w — —1. Hence, if we define /{ = 3{w + 1)/(1 + z) and 
/2 = 3{wACDM + 1) / (1 + z) = (since wacdm = -1), we 
see that for the quintessence dark energy models (with 
w > — 1) /{ > /2 and hence /i > /2 (since initially 
p/po = 1 implying that /i = /2 = 0). Thus we have 

P/PO > PACDM / PACDMO = 1- 

Now since the Hubble function is given by _ff = 
Hq y^fljnoil + z)^ + ^Iqp/po, wc See that for quintessence 
dark energy models H should be such that H > 
Hacdm > C. Hence, quintessence dark energy models 
do not undergo the large degeneracies described in sec- 
tion IHI B 3[ for the range of redshifts covered by the su- 
pernovae observations. A similar reasoning, on the other 
hand, shows that ghost dark energy models can suffer 
from similar large degeneracies, since in that case H can 
intersect the C curve. 



V. IMPORTANT SUBSETS OF 
RECONSTRUCTED CURVES 

In this section we consider different subsets of the re- 
constructed luminosity-distance curves corresponding to 
quintessence dark energy and transient acceleration. 
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A. The quintessence subset 

In this subsection we identify the dark energy with a 
minimally coupled scalar field and reconstruct the so- 
called quintessence solutions defined by w > — 1 (see @ 



Figure 11: Reconstruction of the scalar field kinetic energy for a 
quintessence dark energy model with the corresponding 1, 2 and 3cr 
confidence contours. 
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Figure 12: Reconstruction of the scalar field potential for a 
quintessence dark energy model with the corresponding 1, 2 and 3cr 
confidence contours. 



Fi gure 14: Reconstruction of the EOS with the Universe undergoing 
a transient acceleration with the la confidence contour The upper solid 
line represents the w — —1/3 which demarcates between an accelerated 
and a decelerated expansion. 



B. The transient subset 



VI. CONCLUSIONS 



We also find a subset of reconstructed luminosity- 
distance curves which correspond to models with tran- 
sient accelerating phases - with the universe presently 
decelerating. Figures [13] and [M] show examples of such 
reconstructions for d'/ and the EOS. We have plotted the 
curves with < 118 in order to show that such curves 
fit the data quite well. 

The results also show that to be compatible with the 
observational data, the existence of a transient accelera- 
tion phase, with the universe currently decelerating, re- 
quires that should be very small today. 

Interestingly such transient solutions were not found 
in our recent work Q, where an additional assumption, 
namely dj" < 0, was made. The reason is easy to see. 
Given the smallness of d'/ at present (which follows from 
the present reconstructions), then the constraint d'/' < 
would imply that d" would soon become negative and 
stay negative throughout the past, which is excluded by 
our constraint [B] above. In the absence of the constraint 
d'l" < 0, on the other hand, could have increased, thus 
allowing it to be small today and yet positive in the past, 
resulting in a decelerated phase of expansion for small z. 
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Figure 13: Reconstruction of d'/ when Universe undergoes a transient 
acceleration with the corresponding 1, 2 and 3(T confidence contours. 
The solid line represents di for the ACDM model. 



We have studied the effects of generalizing a re- 
cently proposed model-independent scheme by dropping 
the assumption constraining the third derivative of the 
luminosity-distance. We have seen that this has the con- 
sequence of allowing a wider class of EOS to be explored. 

Using the SNLS supernovae data, we have found 
a number of interesting results reconstructing the 
luminosity-distance curves di, subject to constraints [A] 
and [B] above. 

In particular we have found that although the KCDM 
model is compatible with the data within the 2a signifi- 
cance level, it is not at the la level. 

For models whose EOS becomes ghost somewhere in 
the redshift interval [0, 1], large degeneracies can appear 
in the reconstructions of the EOS due to sudden widen- 
ing of the significance contours, for redshifts larger than 
^ 0.6. We have shown that these degeneracies are a con- 
sequence of the fact that the EOS can in these cases take 
very large values or even diverge, due to the fact that the 
Hubble function can approach or intersect the curve C . 
Our results demonstrate that to decrease these degen- 
eracies, the error bars need to be reduced dramatically, 
particularly for data at redshifts larger than 0.6. We 
have made a concrete study of how such degeneracies can 
be reduced and found that constraints on the accuracy 
of data increase with increasing redshifts. Studies of this 
type could allow limits to be put on the extent to which 
the degeneracies on EOS may be reduced in near future. 
Interestingly, we find that this divergence is absent for 
quintessence dark energy models (with w > —1), with the 
consequence that quintessence dark energy models are 
not affected by such degeneracies. In this case ACDM 
model appears to be among the best fitting models. 

We should add that the interpretation of such possi- 
ble blows ups in the EOS will depend on the theoretical 
framework chosen to interpret the observations. Clearly 
such blow ups can occur if the EOS is treated as an ef- 
fective mathematical representation of a theory such as 
the F{R) theories [13 • On the other hand if the EOS 
is treated as the physical representation of dark energy, 
then according to standard physics the blow up needs 
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to be avoided if it corresponds to a blow up in pressure. 
Things are less clear if it corresponds to the vanishing of 
the density at a finite redshift. 

We have also found examples of reconstructed solutions 
showing transient acceleration while fitting the data very 
well. All such curves, however, require a very small to- 
day, corresponding to a very recent phase of deceleration, 
which thus makes them rather artificial. 

These new results highlight the range of new possibili- 
ties that the generalisation to the reconstruction scheme 
considered in this paper allows. 
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